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Abstract - We estimate the cosmological variation of the proton-to-electron mass ratio 
g = m p /m e by measuring the wavelengths of molecular hydrogen transitions in the early 
universe. The analysis is performed using high spectral resolution observations (FWHM m 
7 km/s) of two damped Lyman-a systems at z a b s = 2.3377 and 3.0249 observed along the 
lines of sight to the quasars Q 1232+082 and Q 0347—382 respectively. 

The most conservative result of the analysis is a possible variation of g over the last 
~10 Gyrs, with an amplitude 

Ag/g= (5.7 ±3.8) x lO” 5 . 

The result is significant at the 1.5a level only and should be confirmed by further observa¬ 
tions. This is the most stringent estimate of a possible cosmological variation of // obtained 
up to now. 
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Introduction 


Contemporary theories of fundamental interactions (SUSY GUT, Superstrings/M- 
theory and others) predict that fundamental physical constants change in the course of 
the Universe evolution. First of all, coupling constants vary with increasing energy trans¬ 
fer in particle interactions (corresponding to the so-called “running constants”). This 
effect has been proved by high-energy experiments in accelerators. For example, the fine- 
structure constant a = e 2 /ftc equals 1/137.036 at low energies (E —► 0), but increases 
up-to 1/128.896 for energy E ~ 90 GeV (Vysotsky et al. 1996). Such “running” of the 
constants has to be taken into account in any consideration of the very early Universe. 
Another prediction of the current theories is that the low-energy limits of the constants 
can vary during the cosmological evolution and take different values at different points of 
the space-time. There are several reasons for such variations. In multidimensional theories 
(Kaluza-Klein, “p-brane” models and others) variations of fundamental physical constants 
are a direct result of the cosmological evolution of extra-dimensional sub-spaces. In some 
theories (e.g. Superstrings) variations of the constants are a consequence of the evolution 
of the vacuum state (a vacuum condensate of some scalar field or “Quintessence”). 

Clearly, experimental detection of such variations of the constants would be a great 
step forward in our description of Nature. Recently, Webb et al. (2001) announced the 
detection of a possible variation of the fine-structure constant, A a/a = (—0.72 ± 0.18) x 
10” 5 , over a redshift range 0.5 < z < 3.5. The method used by these authors is based on 
the measurement of the variation of a large number of transitions from different species. 
This decreases significantly the statistical errors. However, the estimate of the systematics 
becomes more complicated than in the method used earlier where only one species was 
considered (e.g. Ivanchik et ah, 1999). In any case, this exciting result should be checked 
using some other method. 


Proton-to-Electron Mass Ratio 

Here we test a possible cosmological variation of A/i//x, where /i is the proton- 
to-clectron mass ratio at the epoch z = 2.3 — 3.0. It should be noted that a variation 
of a, in principle, implies a variation of /j, because any kind of interaction inherent to 
the particle gives a contribution to its observed mass. This means that any variation 
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of the interaction parameters has to produce some variation of the particle mass, and 
consequently g. Unfortunately, the physical mechanism generating the masses of the proton 
and the electron is unknown up to date. Therefore, the exact functional dependence of 
g(a) is unknown too. Nevertheless, there are some models which permit to estimate the 
electromagnetic contribution to the mass of proton and electron (e.g. Gasser & Leutwyler, 
1982) and dependence of m p , m e , and a on a scalar field which may changed during the 
evolution (e.g. Damour & Polyakov, 1994). There are model relations between cosmological 
variation of a and m p (e.g. Calmet & Fritzsch, 2001). In addition, a curious numerical 
relation may be mentioned: the dimensionless constant g = m p /m e approximately equals 
the ratio of the strong interaction constant g 2 / (he) ~ 14 to the electromagnetic interaction 
constant a = e 2 /he ~ 1/137, where g is the effective coupling constant calculated from 
the amplitude of the 7T-meson-nucleon scattering at low energy. 

At present, the proton-to-clectron mass ratio is measured within the relative accu¬ 
racy of 2 x 10 -9 and equals g = 1836.1526670(39) (Mohr & Tailor, 2000). Laboratory 
metrological measurements rule out considerable variation of g on a short time scale, but 
do not exclude changes over the cosmological time, ~ 10 10 years. Moreover, one can not 
reject a priori the possibility that g (as well as other constants) takes different values in 
widely separated regions of the Universe. This should be directly tested by means of astro- 
physical observations of distant extragalactic objects. Measurements of the wavelengths of 
absorption lines in high-redshift quasar spectra is a powerful tool to check directly the pos¬ 
sible variation of the constants over the cosmological time (e.g. /j and a) between present 
days and the epoch at which the absorption-spectrum has been produced, i.e.~ 10 — 13 Gyr 
ago. 

Up to now the most stringent estimate of possible cosmological variation of g was 
obtained by Potekhin et al. (1998), viz Ag/g < (—10 ± 8) x 10 -5 . 

Sensitivity Coefficients 

The method used here has been suggested by Varshalovich & Levshakov (1993) 
and is based on the fact that wavelengths of electron-vibro-rotational lines depend on the 
reduced mass of the molecule. It is essential that this dependence differs for different 
transitions. This makes it possible to distinguish the cosmological redshift of a line from 
the shift caused by a possible variation of g. The change in wavelength A* due to variation 
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of 11 may be described (in the case A/x//x << 1) by the sensitivity coefficient K t defined as 


_ p dA* 
“ Xi d/i 


( 1 ) 


Such coefficients were calculated for the Lyman and Werner bands of molecular hydrogen 
by Varshalovich & Levshakov (1993), and Varshalovich & Potekhin (1995). Thus, the 
measured wavelength A, of a line formed in the absorption system at redshift z a b s can be 
written as 

A i = X° ■ (1 + z a bs ) • (1 + Ki ■ Afi/p ), (2) 


where A° is the laboratory (vacuum) wavelength of the transition. This formula can be 
written in term of redshift Zi — A*/A° — 1 as 


+ Zabs “t“ b " A /' j 


( 3 ) 


where b = (l + z a bs)-Ap/p. In reality, z t is measured with some uncertainty which is caused 
by statistical errors of the astronomical measurements A i and also errors of the laboratory 
measurements of the wavelengths Xf Thus, a linear regression analysis yields z a bs and b 
(consequently A p/p) with their statistical significance. 


Observations and Results of Analysis 

To measure the possible variation of p we use high resolution spectra (FWHM « 
7 krn/s) of quasars obtained with LIVES at the 8.2-m ESO VLT Kueyen telescope. Two 
H 2 absorption systems were analysed at z a b s = 2.3377 in the spectrum of Q 1232+082 
(Petitjean et ah 2000), and at z a b s = 3.0249 in the spectrum of Q 0347-382 (UVES 
commissioning data, see D’Odorico et al. 2001). 

Absorption system of H 2 at z=2.3377 in the spectrum of Q 1232+082 
More than 50 lines of molecular hydrogen (with S/N ~ 10-14) are identified in the range 
3400-3800 A. We have carefully selected lines that are isolated, unsaturated and unblended. 
In this system only 12 lines meet all these conditions (see Table 1). The accuracy on the 
observed wavelengths A * calculated in accordance with Eq. (A14) from Bohlin et al. (1983) 
which takes into account the number of points within the profile of the spectral line, the 
spectral resolution and the S/N ratio. The average uncertainty of the determination of the 
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line centers is a(A,;) ~ 5 mA. For the laboratory wavelengths we have used two independent 
sets of data: A° [M\ from Morton & Dinerstein, 1976; and A° [A] from Abgrall et al., 1993 
(see, also Ronc-in & Launay, 1994). Results of the linear regression analysis of z t -to-K t 
are shown on Fig. 1 for both sets of laboratory wavelengths. They are the following: 
A p/ix = (14.4 ± 11.4) x 10“ 5 [A], and A////X = (13.2 ± 7.4) x 10" 5 [M], 

Absorption system of H 2 at z=3.02f9 in the spectrum of Q 034 7-382 
For the first time, this H 2 system was found and investigated by Levshakov et al., 2001. 
More than 80 lines of molecular hydrogen (with S/N ~ from 10 to 40) can be identified 
in the wavelength range 3600-4600 A. We have reanalysed this spectrum independently. 
For our analysis 18 lines of H 2 were selected which satisfied to the conditions mentioned 
above. The average uncertainty of the determination of the line centers is cr(Aj) 10 mA. 
Parameters of the lines are presented in Table 2. The results of the linear regression 
analysis of Zi-to-Ki are shown on Fig. 2 for both sets of laboratory wavelengths. They are 
the following: A p/p = (5.8±3.4) x 10“ 5 [A], and A////j = (12.2±7.3) x 10“ 5 [M], It should 
be mention that three points on the bottom panel of Fig. 2 depart from the regression line 
more than 3a. Two of them corresponding to L 9-0 R(l) and W 1-0 R(l) transitions have 
the laboratory wavelengths marked by Morton & Dinerstein (1976) as a blended line and 
a line under weak continuum. The third point corresponding to W 3-0 Q(l) transition has 
deviations on the both panels. It may be a result of undetectable blending in the quasar 
spectrum. We do not reject these points because all of them satisfy to the above formulated 
conditions for the line selection from quasar spectra. 


Combined Analysis 

The combined analysis of the H 2 lines from the two systems discussed above allows us 
to increase the statistical significance of the estimate because of increasing the number of 
lines involved and, what is more essential, broadening the interval of K-values. The results 
of linear regression analysis of Q as a function of K t for all 30 lines from both systems are 
shown in Fig. 3. Here Q is the reduced value of the line redshift: 

A _ z i ~ % 

Si 1 1 — ? 

1 ± Z 

where ~z is z a bs corresponding to the absorption system and a particular set of laboratory 
wavelengths. 
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As a result of the combined analysis we obtained the following estimates (for two 
different sets of laboratory wavelengths): 

Ayu//i = (5.7 ± 3.8) x 10 -5 [A ], 

A n/fj. = (12.5 ± 4.5) x 10“ 5 [M\. 

The statistical uncertainties of the laboratory wavelengths are about 1.5 mA corresponding 
to an error of Ap/fi about 2 x 10~ 5 that is in agreement with the errors found from the 
regression. 


Conclusions 

The above results may be considered as a glimpse on possible cosmological variation 
of p. Additional measurements are necessary to ascertain the conclusion. 

In any case, we have obtained the most stringent estimate on a possible cosmological 
variation of p between redshift zero and redshifts 2-3. 

In order to improve the result, it is necessary to measure more H 2 absorption systems 
at high redshift. The most suitable quasars for such analysis are PKS 0528-250, Q 0347- 
382, and Q 1232+082. Their observations with high resolution (FWHM ~ 7 km/s) and 
high S/N ratio (> 30) will strengthen the conclusions. 

In addition, better accuracy of laboratory wavelengths is also desirable because 
the contribution of laboratory statistical errors are comparable to the statistical errors of 
astronomical measurements. 
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Table 1: H 2 lines of absorption system at z a b s = 2.3377 in Q 1232+082 spectrum 


Lines 

A?, A [M] 

A?, A [A] 

\i, A 

A 

Ki 

L 0-0 P(3) 

1115.896 

1115.895 

3724.543 

0.005 

-0.01479 

L 0-0 R(3) 

1112.584 

1112.583 

3713.480 

0.005 

-0.01178 

L 0-0 P(2) 

1112.495 

1112.459 

3713.179 

0.005 

-0.01170 

L 1-0 P(4) 

1104.084 

1104.084 

3685.093 

0.009 

-0.01154 

L 2-0 P(3) 

1084.559 

1084.562 

3619.947 

0.005 

-0.00098 

L 3-0 P(4) 

1074.313 

1074.314 

3585.740 

0.007 

0.00122 

L 3-0 R(4) 

1070.898 

1070.899 

3574.344 

0.006 

0.00439 

L 3-0 P(3) 

1070.142 

1070.138 

3571.834 

0.005 

0.00511 

L 3-0 R(3) 

1067.478 

1067.474 

3562.948 

0.005 

0.00758 

L 3-0 P(2) 

1066.901 

1066.899 

3561.002 

0.005 

0.00812 

L 4-0 P(4) 

1060.580 

1060.580 

3539.908 

0.005 

0.00685 

L 4-0 P(2) 

1053.281 

1053.283 

3515.556 

0.005 

0.01369 






Table 2: H 2 lines of absorption system at z a bs = 3.0249 in Q 0347-382 spectrum 



Lines 

A?, A [M] 

A", A [A] 

A?,, A 

cr(Aj), A 

Ki 

L 

2-0 R(l) 

1077.698 

1077.697 

4337.614 

0.010 

0.00535 

L 

3-0 R(2) 

1064.995 

1064.993 

4286.483 

0.015 

0.00989 

L 

3-0 P(l) 

1064.606 

1064.606 

4284.924 

0.006 

0.01026 

L 

3-0 R(l) 

1063.460 

1063.460 

4280.313 

0.010 

0.01132 

L 

4-0 R(3) 

1053.976 

1053.977 

4242.144 

0.010 

0.01304 

L 

4-0 R(2) 

1051.497 

1051.498 

4232.175 

0.020 

0.01536 

L 

6-0 R(3) 

1028.986 

1028.983 

4141.571 

0.015 

0.02262 

L 

7-0 R(l) 

1013.434 

1013.436 

4078.977 

0.007 

0.03062 

W 

0-0 Q(2) 

1010.941 

1010.938 

4068.911 

0.010 

-0.00686 

W 

0-0 R(2) 

1009.030 

1009.023 

4061.215 

0.015 

-0.00503 

L 

9-0 R(l) 

992.022 

992.013 

3992.754 

0.010 

0.03796 

W 

1-0 Q(2) 

987.978 

987.974 

3976.492 

0.010 

0.00394 

W 

1-0 R(l) 

985.651 

985.636 

3967.087 

0.007 

0.00626 

L 

10-0 P(l) 

982.834 

982.834 

3955.814 

0.010 

0.04053 

L 12-0 R(3) 

967.674 

967.675 

3894.798 

0.008 

0.04386 

W 

2-0 Q(2) 

967.278 

967.279 

3893.194 

0.010 

0.01301 

W 

2-0 Q(l) 

966.097 

966.094 

3888.423 

0.007 

0.01423 

W 

3-0 Q(l) 

947.425 

947.422 

3813.255 

0.008 

0.02176 






Q 1232+082 




Figure 1: Results of Zi-to-Ki regression analysis for 1+ lines at z a i )S = 2.3377 
in the Q 1232+082 spectrum. Upper panel: laboratory wavelengths are 
taken from Abgrall et al. (1993). Bottom panel: laboratory wavelengths 
are taken from Morton & Dinerstein (1976). 


















Q 0347-382 




Figure 2: Results of Zi-to-Ki regression analysis for H 2 lines at z a i )S = 3.0249 
in the Q 0347-382 spectrum. Upper panel: laboratory wavelengths are taken 
from Abgrall et al. (1993). Bottom panel: laboratory wavelengths are taken 
from Morton & Dinerstein (1976). 









































Figure 3: Combined regression analysis of Q-to-K t for the systems at 

= 2.3377 (triangles) and 3.0249 (squares). Upper panel: laboratory 
wavelengths are taken from Abgrall et al. (1993). Bottom panel: labora¬ 
tory wavelengths are taken from Morton & Dinerstein (1976). 



































